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Abstract. We consider the Frohlich ./V-polaron Hamiltonian in the strong coupling limit 
and bound the ground state energy from below. In particular, our lower bound confirms 
that the ground state energy of the Frohlich polaron and the ground state energy of the 
associated Pekar-Tomasevich variational problem are asymptotically equal in the strong 
coupling limit. We generalize the operator approach that was used to prove a similar result 
in the N = 1 case in ILT| and apply a Feynman-Kac formula to obtain the same result for 
an arbitrary particle number N > 1. 



1. Introduction and Main Results 

The polaron was modelled in 1937 by Frohlich in order to explain the phenomenon of 
electrical breakdown in crystals [Ft]- The model describes a single electron interacting with 
the phonon field of the lattice of atoms of a polar crystal. In the time since Frohlich's initial 
studies, the polaron has found interest not only in condensed matter theory, but also as a 
simple example of non-relativistic quantum field theory. 'Polaron' refers to the quasiparticle 
of the electron together with the phonons created by its interaction with the polar crystal. We 
refer the reader to [AD] and [M] for more details. 

The iV-polaron is a generalization of the single electron case, and describes iV electrons 
interacting with the phonon field and also each other through the Coulomb force. There is 
an asymptotic formula for the ground state energy in the strong coupling limit. In the case 
N = 1 it was proposed by Pekar in [P] and proven by Donsker and Varadhan in |DVj . We 
are interested in proving this formula for polarons of an arbitrary electron number, i.e., when 
N>1. 

The Frohlich Hamiltonian describing a single electron interacting with a polar crystal is 
given by 

=p 2 -^(j>(x) + Hf, (1) 

and acts on the Hilbert space L 2 (R 3 ) <g) J 7 , with T the bosonic Fock space over L 2 (R 3 ). Here 
p = — iV is the electron momemtum operator and the phonon field energy is 



H f = / dka*(k)a(k) (2) 

where a*{k)/a{k) are the creation/annihilation operators for a phonon of momentum k. The 
interaction of the crystal modes with the electron is 

1 r r\k 

* ) ^/,.w (t "'* ,+e " a ' |l,) ' <3) 

The coupling constant a > describes the strength of the interaction between the electron and 
the polar crystal. For a careful definition of H^ 1 ' as a self-adjoint operator that is bounded 
from below, we refer the reader to [N] (see also [MS,). For our methods, this definition is not 
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very important as one can always interpret H^ 1 ' as a quadratic form that is bounded from 
below. We define its ground state energy by E^(a), where we have denoted its dependence 
on a. 

The Hamiltonian describing N electrons in a polar crystal (the 'iV-polaron' or 'multipo- 
laron') is given by 

N 

H u ] = E (Pi - + H f + UV C {X), (4) 

<=i 

with X = (x%, xn) £ R 3N . Here p, = — iV Xi is the momentum operator for the i-th electron 
and 

^> = £^, (5) 

is the Coulomb repulsion between the electrons. The Hamiltonian Hjj 1 ' acts on L 2 (R 3N ) ® T . 
Physically, the parameter U is the square of the electron charge and satisfies U > 2a. However, 
we will also consider the regime U < 2a. 

The ground state energy of is denoted by E^\a), that is 

E$\a)= inf {m { u ] \^)- (6) 

11*11=1 

It is well known that E^\a) > — oo and this is even a simple consequence of our methods. 
In fact, when U > 2a, the system is stable; that is, the ground state energy is bounded below 
by a constant times the particle number [FLSTj . 

The strong coupling limit, in which a — > oo, was first studied by Pekar in the 1950s [P]. 
He hypothesized that in this limit, the asymptotic ground state would equal a product state 
ip ® £, for ip an electronic wave function and £ a phonon wave function. If one computes 
(ip C3> ^\H^\ip (g) £) using the Pekar ansatz, it is easy to determine the minimizing phonon 
wave function £ for a given electronic wave function ip. With this choice of £, one is led to the 
iV-particle Pekar- Tomasevich (PT) functional defined for ip £ i7 1 (R 3Ar ) as 

rg°m = £f |v^| 2 dx + C /W M-^dx-a// ^MMyl dxdy , 

TTlJw T^Jr™ F«- x i\ i )r3 x r3 \x-y\ 



i<3 



(J) 



where 



Pjjjx) = ^ / \ip(xi,...,Xi-i,x, x i+ i,...,x N )\ 2 dx 1 ...dx i ...dx N . (8) 

The hat indicates that dxi is omitted in the integration. We define the ground state energy of 
the PT functional to be 



£f>(a) = inf {pff'W : ^ e ^(R 3JV ), \M\ L > = l} • (9) 



In the case iV = 1 we drop the subscript and write £ W(a). This case was studied by Lieb [L] 
who proved that there is a unique minimizer up to translations. The PT functional in the case 
N > 1 is studied in [Lej . where the existence of minimizers is proven. Lewin also addresses 
the problem of binding in the PT functional. In the regime U < 2a the ground state energy 
of the PT functional is studied in |BBj ; they also sketch the derivation of ([7]) from the Pekar 
ansatz. 
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A scaling argument gives 

4^(a)=fig£(l)o». (10) 

Furthermore, because one can construct wave functions where the electrons are arbitrarily far 
apart from one another, the ground state energy of the PT functional is subadditive: 

4^(a)<4%)+4^(a). (11) 

Pekar's hypothesis amounts to the conjecture that in the strong coupling limit, the asymp- 
totic formula 

lim a- 2 E^ N \a)=£i N \l), (12) 

a— >oo 

holds. Here, v = U/a is a fixed constant. 

Because it arises using a variational ansatz, £^\a) is automatically an upper bound for 
Ejf (a). In 1983, Donsker and Varadhan first confirmed (|T2|) in the case N = 1 |DVj . However, 
their proof did not give a rate of convergence and did not easily generalize to other settings. In 
1995, Lieb and Thomas developed a simpler proof using operator methods [LT . They bound 
E^(a) from below by £^(a) minus an error term which is negligible in the strong coupling 
limit, thus proving (IT2"1) . 

In the paper [MS , Miyao and Spohn treated the strong coupling limit in the bipolaron 
case (their proof lies in the appendix) using the methodology of |LT] . It was also claimed in 
[FLST] that the proof in |MSj generalizes to the N > 2 case. Unfortunately, the proof in [MS 
is slightly incomplete. While not confirming the formula (fT2")l . their proof does give a lower 
bound on the energy of wave functions that describe electrons which are all localized to the 
same small box in space (the case that their methods do not address is when the electrons are 
very far apart from one another). Furthermore, this generalizes easily to N > 2. In fact, this 
generalization is a crucial element in our own proof. 

It is our contribution to complete the proof using a Feynman-Kac formula, and to make the 
generalization to the N > 2 case. In particular, we prove a lower bound for the ground state 
energy E^(a) which implies (fT2|) for any N. 

The astute reader will have noticed that we do not impose any symmetry constraints on the 
electron wave functions. Ideally, one would like to obtain a lower bound for the ground state 
energy E^\a) (taken over only antisymmetric wave functions) in terms of £jj f \a) but with 
the infimum in ((9|) replaced by an infimum over the antisymmetric wave functions in iJ^R 3 ^). 
However, our methods do not directly generalize to this setting. 

We now state our main result. In particular, it confirms (fl2"|) . Recall that the ratio v = U /a 
is fixed. 

Theorem 1. (i) For any given v > 2 and e > there exists a constant C t . v > such 
that for allN € IN and a > eN 4 , 

E^ N \a) > a 2 £i N \\) - C^a 9 / 5 N 9 ' 5 . (13) 

(ii) For any given v > and e > 0, there exists a constant D t v > such that for all 
N G N and a > e, 

E ( v N \a) > a 2 £i N \l) - D^a i2 ' 2Z N\ (14) 

The asymptotic formula (fl"2"|) allows us to extract information about the Frohlich Hamilton- 
ian from the PT approximation. In [Lej , it was proven that for any N, there is a v(N) > 2 so 
that the binding inequality 

£l N Hl)<£i N - k \l)+£i k \l), (15) 
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holds for any v < v{N) and k. This clearly implies 

Corollary 1. For any N there is a v(N) > 2 so that for every v < v{N), there is an 
a(N, v) > so that the binding inequality 



holds for any a > a(N, v) and any k. 

In Theorem [I] we have stated two lower bounds which are valid for different relative sizes 
of a and N; both imply (|12[) . It is an eventual goal to prove a lower bound of the kind (when 
v > 2 and hence the system is stable), 



Such a lower bound would imply that one can, in a certain sense, 'commute' the two limits 
N oo and a — > oo when considering the quantity 



In this sense, (fT5)) is not optimal in the N dependence of the error term. The N 9 / 5 term is 
close; however, as stated, the bound does not necessarily hold if one takes N —> oo first. We 
hope that (TT3l) is a first step towards (fT7|) . 

Let us comment on the restriction on a in Theorem l(i). In the physical regime v > 2 
where the PT energy satisfies £u(l) = O(N), the assumption a > eiV 4 is necessary for 
the error in (IT51) to be small. What we mean is that, if one rewrites the RHS of (fT5)) as 
a 2 N[£i N \l)/N — C e ^a _1 / 5 iV 4 / 5 ], then it is easy to see that the error term will be small 
relative to the PT energy only if a ^S> N 4 . Therefore, (fT3f is only useful if the assumption 
a > eiV 4 is satisfied, and so this assumption isn't very restrictive anyway. 

Comparatively, in the regime v < 2, we have that £\f \a) = 0(N 3 ) [BBj . Rewriting the 
RHS of (33]) as TW^^lVTV 3 - C e a" 4 / 23 ], we see that the corrections to the PT energy 
are small regardless of the relative size of N and a (as long as a ^> 1). 

I.I. Outline of Proof. Our proof generalizes the operator methods developed in |LTj . which 
treated the N = 1 case, to polarons of arbitrary electron number. Using a continuous version 
of the IMS localization formula, we localize the electrons to cubes of a fixed side length in R 3 . 
While we have localized each of the electrons to their own cube, we do not, a priori, know 
anything about how far the electrons are from one another. To quantify how spread out they 
are, we partition the electrons into disjoint clusters. Electrons in a single cluster are not too 
far from each other, and the electrons in different clusters are separated by some minimum 
distance. 

The generalization of the methods in [LTj (and also [MS] ) provides a lower bound on the 
energy of the electrons that depends on how far apart the electrons are from one another; 
that is, if they are very spread out, then the lower bound is not very useful. What we would 
like to do is ignore the interaction between electrons lying in different clusters, and apply the 
methodology of |LT] to obtain a bound on the energy of each individual cluster. Since the 
electrons in a single cluster will be tightly packed, the bounds given by |LT] will be good 
enough to handle the energy of the clusters. 

In |FLST| , the Feynman-Kac formula was applied in order to bound the interaction of two 
localized electrons sufficiently far from each other. We generalize this to the A^-polaron case 
to show that we can ignore the interaction between different clusters (at an appropriate error) 



E^(a)<E^- k \a) + E^(a), 



(16) 



E^\a) >a 2 £l N X\)-0{a 9 '\N). 



(17) 



lim 

N—KX) 




(18) 
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and treat each cluster as its own subsystem. Subsequently, we sum up the bounds on each 
cluster energy obtained using |LT) : this will in turn yield Theorem Q] 

Acknowledgements: The research of B.L. was partly supported by NSERC. The research 
of LA. was supported by the German Science Foundation (DFG). The authors wish to thank 
R. Seiringer for many useful discussions and comments on a draft of this paper. The authors 
are grateful to M. Griesemer and D. Wellig for numerous stimulating discussions and access 
to a preprint of their paper [GWj . 

2. Proof of Theorem 1 

Step 1 (Localization). The first step in proving Theorem [1] is to apply a continuous version 
of the IMS localization formula in order to localize each electron into a cube of finite side 
length. 

Fixi?> and for a, G R 3 , i G {l,2,...,n} let C R {a x , a N ) = (-R/2, R/2) 3N +{ ai , a N ) C 
R 3N . Cn(ai, ajv) is the 3iV-cube of side length R centered at (ai, aw). Let Br be the 
set of normalized wave functions for which the support of their electronic part is contained in 
Cn(ai, ...,cin) for some ai G R 3 (i.e. ip G L 2 (Cn(a)) ® T for some a G R 3N ). The following 
proposition is proven in [LTj in the case N = 1, where the IMS localization formula is applied. 
We generalize this to arbitrary N. 



Proposition 1. For R > 0, 

3NjT 



f(«)> inf M<V>-T- (19) 



Proof. Let e > and $ be a normalized wave function such that (^H^l®) < E < ^\a) + e. 
For real-valued ip G C§° (Cr(0)) set ipy(X) — ip(X — Y). A direct calculation leads to 

f dY{^\H^ N) \^)= I dY\^(Y)\ 2 ((m { u N} m) + \m\U^)- m 

R 3JV J R 3N V / 

We can take ip so that ||V^||^2 < 3Nn 2 /R 2 + e and llV'llia = 1 (i.e., it approximates the 
ground state of the Dirichlet Laplacian). For this choice of 



dY 

R 3N 



(VySliT^lW*) " (E^iaj+W^/R 2 +2e){^V*|Vy*> < °- (21) 



There must be a Y G R 3Ar for which the integrand appearing above is non-positive and 
(ijjYo®\'4'Yo®) is nonzero. Therefore, 

^y ( MH ( u N) \i>Y ^)/^Y ( Mi>Y ^) < E { v N) (a) + 3Nn 2 /R 2 + 2e. (22) 

Since i/)y <i>((V'Y < I ) |V'Yb^ ) ))~ 1 ^ 2 e tne claim is proven. □ 

We now fix a $ G Br for an R > which will be chosen later. Each of the N electrons 
described by <& is located in a cube Qi of side length 7? centered at some a, G R 3 . The remainder 
of the proof constitutes finding a lower bound for ($|.ff£^|<I>) which holds independently of 

i.e., finding a lower bound for the RHS of ifT^j) . 

Step 2 (Partitioning of electrons into clusters) . The direct generalization of [LTj and [MS] 
to the 7V-polaron case produces a lower bound on (<i>|.ff ( y(a!)|<i>) which depends on how far 
apart the localized electrons are from one another. Herein lies the main difficulty in our proof: 
we do not, a priori, know that the cubes {Qi}fL 1 within which the electrons are contained are 
close enough to each other so that this bound is useful. In order to overcome this challenge, 
we must quantify how spread out the electrons are. This is accomplished by partitioning the 
electrons into clusters. 
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Let d(i,j) = dist{Qi, Qj} be the distance between the cubes containing the supports of the 
i-th and j-th electrons. Let d > be a fixed constant. We define the relation ~ on {1, N} by 
j k iff 3 a sequence %m) with %\ = j, %m = k satisfying d(ii, h+i) < d for 1 < I < M. 

Clearly ~ is a relation of equivalence and we let G±, Gi be the disjoint equivalence classes 
to which we will give the suggestive name of 'clusters.' The interpretation of ~ is that we 
have placed two electrons in the same cluster iff either they are less than d apart or they are 
connected by a sequence in which consecutive electrons are less than d apart. 

There are two important properties of these clusters: 

(PI): Let \Gi\ — Ni be the cardinality of the cluster Gi. Then it is clear that the union of the 
supports of the electrons in Gi is contained in a cube in K 3 of side length no greater 
than Ni(R + d). 

(P2): If the electrons i and j are in two different clusters, then d(i,j) > d. 

Property (P2) will allow us to apply the Feynman-Kac formula to bound the interaction 
between electrons in different clusters. Property (PI) will allow us to apply the methods of 
|LT| to bound the remaining energy of each cluster of electrons. 

Step 3 (Removing the inter-cluster interaction). For the cluster Gi, we define 

E(G i )=in£ , (<p\HV f,) \<p), (23) 
<p 

where the infimum is taken over Ni— polaron wave functions <p which have the same support 
properties as the electrons in Gi. More precisely if the electrons in Gi are localized in the 
cubes Qj 1 , Qj N . then the <p appearing in the above infimum are also supported in the same 
Ni cubes (i.e., ip 6 L 2 (Qj 1 x ... x Qj N .) <8> F). Physically, E(Gi) may be interpreted as the 
lowest energy the electrons in the cluster Gi can have (ignoring, of course, the existence of the 
other N — Ni electrons) . 

The following lemma provides a lower bound for Ey (a) in terms of the cluster energies 
{E(Gi)}i and an error term due to the inter-cluster interaction. 

Lemma 1. For our fixed <& € Br we have, 

($, # <w#> > y E(Gi) y y ( u ^ ) . (24) 



kiy£k 2 i£G kl 



In particular, we have, 



<*,FW*)>£%7,)-^. (25) 

i 

Ifv>2 andd> A^/ZR/(v - 2) then, 

WnSPw >£%?<)• (26) 

i 

Remark 1 . This is a direct generalization of Lemma 1 in [FLSTj , in which the N — 2 case is 
proven. 

Proof. The Feynman-Kac formula [R] implies that the infimum of the LHS of over nor- 
malized wave functions localized in Q\ x ... x Qn equals 

- Urn ^logZ Ql ,..., Qjv (T), (27) 

1 — >OQ 1 
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where 



AT) 



x exp 



da- 



dxi... / dx N / dW Xl (u 1 )..AW XN (u N )xQ 1 (ui)---XQ N (uN) 
i JQn J 

r * uy( T * ' 



i<3 



(28) 



Above, dWj. denotes the Wiener measure of closed Brownian paths : [0, T] —> R 3 , starting 
and ending at . XQj ( w j ) is a characteristic function which is 1 if uij stays inside the cube 
Qj for all times and otherwise. 

We use our choice of clusters to bound the integrand in the above equation. Specifically, 
we obtain for paths lying in the cubes {Qi}i (the only paths for which the characteristic 

functions XQ«( w i) are nonzero), 



-Ul N 

d A^E 

2 ^ 

*i3= 



elf 



»<3 



df 



< 



£ 

*!=i 



a I ds- 



E 



dt 



+ E^E 



o - Wj(i + s) 

f7 



o K(t)-w,-(t)| 
, Jo 



dt 



i<j 



Wi(t) - Wj(t)\ 



J 



d(i,j) 2d{iJ) + 4V3R 



(29) 



The terms appearing in the first line of (|29p involve electrons either lying in the same cluster or 
lying in different clusters; in the subsequent upper bound, the former have been left unchanged, 
but the latter have been bounded above using what we know about the inter-particle distance. 
To be precise, we have used that d(i,j) < \uii{t) — ojj(t + s)\ < d(i,j) + 2\/3R for all t and s. 
It follows that 



1 1 1 

-logZ Qu ... tQN (T)<J2^ogZ Gk (T)- £ £ 



fe=i 



jec k2 



U 



d(i,j) 2d(i,j)+4V3R 



(30) 



where, if Gu = {ii, ...,ijv fc }, then 



Z Gk {T) := Z Qtl ,..., Q (T), 



(31) 



and ZQ ii> ... t Q i (T) is defined in the obvious way by (|2"5|) . Taking the limit T — > oo on both 
sides of (001) yields {H}. Since J7 > 0, (J2SJ follows from (JSJl. Solving for the d that makes 
the terms appearing on the RHS of (|2"4"| vanish and recalling property (P2) of the clusters we 
obtain (J2BJ). □ 



Step 4 (Lower bound for cluster energies). Lemma[T]tells us that we need only bound each 
E(Gi) from below. To obtain such a bound, we will follow very closely |LT] and [MS , and 
our methodology is essentially a straightforward generalization to the N > 2 case. We give 
here the proof for reader convenience, and also refer those interested in more details to these 
papers. 
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Our first step in deriving the lower bound will be to replace the Hamiltonian with 
an ultraviolet cutoff Hamiltonian. For Ki > (the value of our ultraviolet cutoff) let 



H 



N it U 
Ki 



1 



8N ia 



(e tkx *a(k) + e~ tkx *a*{k)) 



U 



dka* (fc)a(fc) + > J \ 
\k\<Ki \ x i - x j\ 



(32) 



We now prove: 



Lemma 2. Let be the Ni-polaron Hamiltonian as in and H^ u be the cutoff Hamil- 

tonian as in 132\). Then, 

T N it U _ 1 
2' 

Remark 2. By completing the square, the phonon variables can be eliminated from H^' U as in 
[LTj . for some cut-off K. This leaves only the positive kinetic energy, positive Coulomb energy 
and a negative constant proportional to K. Therefore, applying Lemma [2] to shows that 

E ( v\a) > -co. 

Proof. For notational simplicity we write N and K instead of Ni and Ki, and relabel the 
electrons in the cluster G; as 1, ...,N. We define for each electron the vector operator = 
(Z» where 



13 17 - "jfj 



(33) 



fc|>K 



dfc y fc,ai , 



(34) 



A standard calculation yields 

N 3 



t=l J=l 



AT 



dfc 



jc V3»r|*| 



e* x *a(fc) + e-^'^o'ffe)] =: flj, (35) 



where the rightmost equality defines Hj. Above, Pij is the i-th electron momentum in direction 
j. We now bound the LHS of ([331) . Fix e > 0. For arbitrary normalized ip, 



j^{ V \[ Pij ,zf -Zf>] \<p) 

3=1 



< 2(^ 2 |^) 1 /2 (¥ ,| _ (Z (o _ z^fW 2 

< 2(tp\p 2 l \(p) 1 / 2 (ip\2(Z^*Z^ + ZWZW*)^) 1 

< e(ip\rf\(p) + ^(ip\Z^*Z^ + Z^Z^*\<p) 
= e{<p\p 2 W) + -{ V \Z^*Z^\ V ) ■ ' 4(> 



/2 



In the last line we have used 



3=1 



7 {i)* 
Z J ' Z 3 



Additionally, the standard upper bound 



\k\>K 



2a 



dfc 

27T 2 |fc| 4 



2a 
~txK 



Z (i)* z d) < ^7V >if , 



holds, with 



dfca*(fc)a(fc). 



(36) 
(37) 

(38) 
(39) 



\k\>K 
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At this point, we take e = 8Na/irK. With Hi defined as in (|35]). the bounds (|36|) and (|38| 
imply 

MHjM < + (<p\*>k\<p) + 5- (4°) 

from which it follows that 

< } - = ^E^+ N >« - ^ * 4 (41) 

The above inequality proves Lemma H □ 



Let us focus on a single cluster G{. We prove the following bound for the cluster energies: 

Proposition 2. For the cluster Gi let Pi and Si be positive constants, with 2Si < 1. Assume 
that the ultraviolet cut-off Ki satisfies 8aNi < irKi. Then, 

E(Gi) > Fm : = o^(l) 1 8jV - g^±jM^ _ + 1} a _ I 

V V »; ^ u (!_ 2^)2(1 _ ^t) 7r<5< l P ^ 2 

(42) 

with Vi = U(l — 25i)/a. 

Remark 3. The meaning of the constants P, and Si will be made clear below. The values of 
all the parameters introduced in Proposition [5] will be chosen later (each claim of Theorem [T] 
requires a different choice of constants). Our definition of P(-/V,-) is suggestive; we are thinking 
of the lower bound for E(Gi) as a function of the particle number N t . Indeed, we will later 
see that the parameters 5$, Ki, etc. will depend only on N t (and a, of course). 

Proof. For notational simplicity of the proof we will omit the subscripts where the meaning 
is clear and write N instead of Ni, K instead of Ki, etc. In order to get a lower bound on 
E(Gi) we fix an arbitrary normalized wave function with the required support properties. 
By Lemma [2] it suffices to consider the cut-off Hamiltonian as defined in (|3"2")) . We will find 
a lower bound for (\E , |i/^' C/ |\E') which is independent of ^ . Recall that our ultraviolet cut-off 
is denoted by K > 0. We will now split the sphere Bk '■= {k £ R 3 : |fc| < K} into cubes of 
length Pi = P (whose value will be determined later). For n = (m, 712,113) £ Z 3 we define 
D P (n) = [mP - P/2, mP + P/2] x [n 2 P - P/2, n 2 P + P/2] x [n 3 P - P/2, n 3 P + P/2] C R 3 , 
the cube with center Pn and side length P. Let Ap = {neZ 3 : Dp(n) n Bk 7^ 0}- For the 
cardinality of Ap, denoted |Ap|, we have 

9K 

|Ap|<(— + 1) 3 . (43) 

For n £ Ap set B(n) — Dp(n)DBK , and let fcs(n) be any point in B(n). The electrons described 
by VP are localized in a cube of side length N(R + d). By applying a unitary transformation 
(i.e., translating the system) we may assume that the cube is centered at the origin. Then, 
using that the diameter of a cube is \/2> times its side length, we obtain for Xi in the support 
of * that 



3P/V 

\jk. Xi _ e ik B(nyXi] < | fc _ fcs(n) | N < <l_ (R + d y (44) 
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For S — Si > 0, we have by completing the square, 

JV 



neA r 



E -v^E 



dk 



~~[JB(n) V2ir\k 
dka*(k)a(k) 



(e ik - Xi - e ikB ^r Xi ) a(k) + (t 



> 



B(n) 

aN 
2^5 



E 



\kV^ ' 



1=1 



2n5 



Therefore, 



*|^' U |*)>(*|^({fc s( „)})|*}- 



9aiV 4 P 2 (i? + d) 2 A: 
2^ ' 



(45) 



(46) 



with 

H K ({k B (n)}) = (1 



8aiV,A 2 rrV ^ 1 



irK 



i=l 



\Xi Xj\ 



E 



(1-6) 



N 

[ dka*(k)a(k)- V^V / / fc (e' fcB <")- a: 'a(fc) + e- ik w Xi a*(k)) 

JB(n) i=1 JB(n) V2n\k\ 



Define the block annihilation operator A n for n G Ap by 

-1/2 

/ /' (I/.- \ 

A-n = 



dk 

B(n) 2^\kJ 



dk 



B(> 



V2w\k 



■a(k). 



(47) 



(48) 



Each A n is a normalized boson mode satisfying [j4„,A*/] = 5 nn '. By the Schwarz inequality, 
A* n A n < J B{n) dka*(k)a{k). Therefore H K ({k B ( n )}) > H§ 1 ({k B{n) }), with the block Hamil- 
tonian defined by 

A 



af ({k B(n) }) = (l 



8aN 
ttK 



)E^ + c/ Et^^ 

trf ~f-\ x i- x j\ 



E 

n£Ap 



(i - 



Bin) 2^ 2 |fc| 2 



1/2 jv 

j (e ik "or Xi A n + e- ikn ^r^A' n ) 



In summary, we have proven 



9aN i P 2 (R + d) 2 K 
2^6 ' 



(49) 



(50) 



Note that the above inequality holds for any choice of {k B ^}. Our goal is now to find a lower 
bound for (^>\H^ 1 ({k B ^})\^) . A lower bound that is independent of ^ will in turn imply a 
lower bound for E(Gi), by the above inequality and LemmaEJ 

We will eventually use coherent states to relate the block Hamiltonian to the PT functional. 
Before doing so we would like to make one observation regarding the block Hamiltonian. 
^K l ({^B(n)}) contains terms involving creation and annihilation operators of only finitely 
many functions. If we denote the span of these functions by .M, then we note that the bosonic 
Fock space T = J"(L 2 (IR 3 )) (where denotes the bosonic Fock space over the Hilbcrt 

space H) is unitarily equivalent to F(M) ® J-(M ± ). Here, we recall that J-{H\ © H 2 ) = 
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Furthermore, because the block Hamiltonian contains no terms associated with the functions 
in ft/t , we can identify it with the operator H^ l ({kB( n )}) ® ^-j r (M ± ) which acts on the Hilbcrt 
space L 2 (R 3N )(g>J 7 (M)(g>J : '{M ± ). Here, we are thinking of H§ l ({k B ( n )}) defined by (@9| as an 
operator acting on L 2 (R 3N ) (g> F(M). Since the operators A®t and A have the same ground 
state energy, it is therefore sufficient to obtain a lower bound for the ground state energy of 
H K l ({ k B(n)}) acting on L 2 (R 3N ) (g> T(M), for a choice of {k B{n) } = {fc(l)„} which we now 
fix. Moreover, since the unitary transformation taking L 2 (R 3N ) ® T to L 2 (R 3N ) ® <S> 
.F(.M ) acts only on the phonon variables, we can still assume that the wave functions that 
ffjf ({fc(l) n }) act on have the same support properties as before - that is, they are localized 
in a box at the origin. 

We now fix an arbitrary wave function in L 2 (R 3N ) which we continue to denote by 

W which has the required support properties. We claim that instead of finding a lower bound 
for the energy of \& using the {k(l) n } we fixed previously, we can instead take the supremum 
of (^\H§- l ({kB(n)})W) over au choices {k B (n)}> at the cost of replacing (1 — 5) by (1 — 26) 
in the definition of ({kgin)}) an d incurring an error term which is equal to the constant 
appearing in (|50|) . What we mean by this is the following. Let {&s( n )} be another choice of 
the fixed points in each B(n). The arguments of (|44|) - (|46|) . applied to the block oscillators 
A#(k) instead of the a#(k), show that 

'({*(l)n})|*) > (mx({kB(n)})\*) S - 9aNip2 ^ + d ^ K . 

(51) 



We can even take the supremum on the RHS of (|51jl over choices of {fc_B(„)}. Define H^ l ({kB( n )}) 
by the RHS of (|49|) except with the coefficient (1 — 25) in front of the block number operator 
instead of (1 — 8). We will now prove the bound 



a 2 £L N) (l) 2K 
sup <*|2tf'({* fl(n) })|*) > - 1 ' - (-p + l) a (Wi 



' T Bl 



Let us see that such a lower bound implies Proposition [5J Since the RHS of (|B"2| docs not 
depend on this, together with (|5ip. implies a lower bound for the ground state energy 
of Hk ({fc(l)„}). Recall that this operator is acting only on wave functions in the Hilbert 
space L 2 (R) (g> that have the aforementioned support properties. As we have previously 

argued, this implies a lower bound for E(Gi). Collecting (|33l) . ([5u]) . (f5"Tj) and (|5^|) we see that 
the bound we obtain is in fact 

E(Gi) > F(N), (53) 

with F(N) defined as in (H^)) . which proves Proposition [51 

We have left only to prove (|55|). For a given £ = (£ n )neA P € C' Ap we introduce the coherent 
state in F(M) defined by 

|0 := [] et^-e^SlM, (54) 

neAp 

with fi^f the vacuum state in T(Ai). With this definition, we have that 

/ dedriexei = / wucxci = / dcdr(ie«i 2 - i)io<ei = kk, (55) 

where d£d£* = rineAp 7I " -1( ^nd^n anc ^ denotes Lebesgue measure. 
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Let := (&,£)f(m) S L 2 (R 3N ) be the wave function obtained by integrating out the 
phonon variables. With the identities (|55]l we have, 



(56) 



where 



h ({**(«)}) = (1 



8aN s 
irK ' 



e^+^e,^-^ 

<—* ^— ' \Xi - Xj\ 



E 

ngAp 



(l-2J)(|e„| 2 -l)-^ 



dk 

2n 2 \k\ 2 



1/2 

E(e» 



i=l 



Completing the square we obtain 



JV 

<*el* ({*«(»)}) I*?) > <*el(i - ^) Ep? + ^E 



t<3 



E 

n€Ap 



47ra 
1-25 



\k( k B(n))\ 2 
B{n) \k\ 2 ||*f|l^ (R 3JV) 



dfe- 



(1-2,5)11^1 



L 2 (R 



3ivi |Ap 



where p^{k) = (2ir) 3 / 2 J" R3 dxe lkx p^{x) is the Fourier transform of 

N r ~ 
Pd x ) = E / |*{(a;i, ...,a;i_i,x,x i+ i, ...,a; A f)| 2 da;i...dx l ...dxAf. 



We have the inequality 



inf 



dfc lft(**w )|a < 



{feB ( „)}7B(„) |fc| 2 ||* e | 



L 2 (R 3JV 



dfc- 

fl(n) |fc| 2 ||*{| 



L 2 (R 3 



and so, 



sup /dedr(*d^ ({*b(u)» l*€> 



L 2 (R 3JV ) 



> 



(*el(i 



47ra 



dA- 



w 1 



(1-25)11* 



) 



Jll_L 2 (R3«) 



ia p | d^dr 



dedni^|l^ (K 3 N) (l-^) 
a 1 



A' 



<*elEi 



i=l 



8aA 



(i-2^)(i-^)||*dlV« 

for ^ = ^/ ||*^ IIl2( R 3jv) • Above, we have used 

Rk)g(k) l 



dxdy 

'xR3 p-y| 



Er^b^> 

-(1-2<5)|A F | 



dfc- 



\k\< 



4tt 



, , f{x)g(y) 

dxdy- r . 

R3 X[R 3 \x-y\ 



(57) 



(58) 



(59) 



(60) 



(61) 



(62) 
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The integrand enclosed in square brackets in the last two lines of (|6Tj) is just a constant 
times the PT functional for <J>£ with the coefficient a replaced by a(l — ^-^-) -1 (l — 2S)~ 1 and 
U replaced by [7(1 — ^y-) -1 • Using this and applying (JTUJ) we obtain the lower bound, 

~ N U 1 

i=l 1 _ TIT i<j l Xl ^1 



2 



^)(l-2fl||* c ||* / J R3xP / Xdy \x-y\ - £ ~» (1) (i-^f) 2 (l- 2 A) 2 ^ 



(63) 



tit >y - * u > n*eiii3(iF 

for z7 = f/(l - 25)/a. The equality (J5B) together with the bounds ([61]) and J63J prove ([52]) , 
thus completing the proof of Proposition [2j □ 
Step 5 (Choice of parameters) Let us summarize the proof so far. In Proposition [TJ we 
proved that the ground state energy Ey(a) is bounded below by the infimum of (^\H^ |$) 
taken over wave functions in Br (the set of normalized wave functions with electrons localized 
in boxes somewhere in space) minus an error term. Subsequently, for an arbitrary wave 
function $ in Br we separated the electrons into the clusters {Gi}i and have shown in Lemma 
[JJthat the energy of $ can be bounded below by a sum of cluster energies, '}2 i E(Gi), minus 
an error term. In Proposition [2] we derived a lower bound on each E(Gi). It remains only to 
collect all of these bounds and choose values for the parameters (e.g., the Ki, R, etc.) that we 
have introduced in the proof. We will subsequently obtain a lower bound for the ground state 
energy. 

In the case U > 2a we take d = y/3R/(i> — 2) and using equations (|19p . (|2"6"|) and (|4"2")l we 
obtain that (recalling the definition of F(Ni)), 

E^(a)>mmJ2(F( Ni )- 3 -^-), (64) 



1=1 



where the minimum is taken over all the possible choices of positive integers N±, Ni, with 
Ni = N, for any I. By taking the minimum, we have removed all reference to the original 
wave function $, and therefore obtain a lower bound for the ground state energy. Recall that 
we are thinking of F(Ni) as a function of the particle number Ni, so that the constants on the 
RHS of (|4*2"j) are thought of as depending on how many particles are in the cluster. 

We will now make our choice of the parameters Ki,R, etc., so that (|64|) yields the claim 
(|13p of Theorem [TJ Let us first explain the reasoning behind our first choice of parameters. 
In the U > 2a regime, the Pekar-Tomasevich energy £jj f \a) is of order N [FLSTj . As a 
consequence, the multiplicative errors appearing in (|42|) (i.e., the 1 — 26 and 1 — ^fjf terms) 
result in corrections to the leading order term a 2 Ei N \l) which are of the order of a 2 NS and 
a 2 N^f. Therefore, due to (02]), the corrections contained in F(Ni) - to the leading 

PT energy term are (up to constants) , 

a 3 N? 2 Ni aNfPfR^K, K? 

aTNiSi, *—! , — §-. (65) 



In order to optimize our bounds, we choose the parameters so that the powers of Ni and 
a appearing in the above terms are all equal. A simple calculation shows that we should 
take Ki = Nf /5 a 6 / 5 c 1 , R = c 2 a- 9 / 10 N' 2 ^ , Pi = c 3 a 3 / 5 Nr 2/5 and Si = c 4 a- 1 / 5 iV i i/5 where 
ci,C2,C3,C4 > are constants. In order to ensure that the hypotheses on Si and Ki stated 
in Proposition [2J are satisfied, we choose the constants to satisfy a > N 4 x max{(2^C4/ \v — 
2)) 5 , (8/(7rci)) 5 }. There is a subtle point here. The leading order term in our lower bound is 
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a 2 £~ N ^ (1), with Pi = (1 — 2<5j)z/, which is slightly smaller than the original v. The assumption 
that a > N 4 (2vci/(iy - 2)) 5 ensures that P t > 2 so that £^ Ni \l) is order A^. 
With this choice of constants, it follows that 

Ml n 9/5 N 9/5 ^-rvo-^ 

■\2 



RHS of m > 5 ff» '<') — , - JW*L**=*£«te 



- X 3 _ a ^ N ^ _/3 6o^\ 

V c 3 J 4 \2 c 3 J 

with v = (\ — 2c4 : a~ 1 / 5 N 4 / 5 )i>. Here we have used subaddivity of the PT ground state energy 
and also the fact that the PT ground state energy is negative. We have also used the elementary 
inequality J^. Nf < N p for p > 1. Note that because it is a lower bound for the ground state 
energy, (p6| alone confirms the asymptotic formula (fT2")h 

We are left with showing that the RHS of (| 13[) is a lower bound for the RHS of Q66p. 
Clearly, is of the same form as (fT5]l . except for the multiplicative errors in the leading 
order PT energy term. The terms appearing on the denominator of the leading order term 
obviously result in corrections that are order iV 9 / 5 a 9 / 5 , leaving us with a lower bound of form 
a 2 £~ N \l) — CW 9 / 5 o; 9 / 5 . The PT energy £~ N \l) (with a weaker Coulomb repulsion) can be 
bounded below by £u (1) (with the original stronger Coulomb repulsion) minus an error term 
which is order (y — v)£i N \l) using the concavity of the function v — > ££^(1)0 This proves 
the first claim of Theorem [TJ 

For general v > the Pekar-Tomasevich energy S^j (a) is always at most of the order iV 3 
BB . To derive (|T3)) . we repeat the above argument, but use ([25|l instead of (|26j) to bound 
Ey(a) by the cluster energies {E(Gi)}i. We choose for our parameters Ki — Nia 27 ^ 23 ci, 
R = c 2 a~ 19 / 23 N~ 1 , d = C5-R, Pi = 03a 13 / 23 and 5, = c^oT 4 ! 23 . With these choices, we obtain 
an inequality similar to (l66l) except for the fact that A r9 / 5 's have become A^ 3 's and that there 
is an additional term arising from our application of the Feynman-Kac formula (i.e., the extra 
term on the RHS of ([25]) ) which contains a higher power of a, namely a 42 / 23 . The same 
argument as that above then yields the RHS of (fT4l) as a lower bound for E^\a). Note that 
we require a > max{(2c4) 23 / 4 , (8/(7rci)) 23 / 4 } in order to satisfy the hypotheses of Proposition 
[2] This completes the proof of Theorem [T] 
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